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The relationship is established between the Berry phase and spin crossover in condensed matter physics 
induced by high pressure. It is shown that the geometric phase has topological origin and can be considered 
as the order parameter for such transition. 
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Spin crossover in condensed matter physics is a 
transformation of a system with one spin Si at each 
lattice site into another state with spin S2 induced by 
some external parameter like strong magnetic field, high 
pressure etc. It accompanies by the energy level E\ 
and E2 crossing, where E a is the local energy of the 
magnetic ion with spin S a (a = 1,2). Recently spin 
crossovers in magnetic oxides have been found under 
high pressure in FeB0 3 [1], CdFe 3 [2J, BiFe0 3 [3], 
Fe 3 04 g]. Below the Curie temperature of magnetic 
order spin crossover is accompanied by the sharp change 
of the magnetization, nevertheless it may be observed 
in the paramagnetic state like in CdFe 3 (B0 3 )i [2] as 
the sharp change of the XES satellite/main peak inten- 
sity ratio with pressure increase. There is no thermody- 
namic order parameter that can distinguish one phase 
versus the other. In this paper we discuss the low tem- 
perature limit and claim that at T = spin crossover 
is a quantum phase transition. The order parameter 
for such transition has topological origin and we calcu- 
late the geometrical phase that characterizes the spin 
crossover. 

Quantum phase transition (QPT) is characterized 
by qualitative changes of the ground state of many body 
system and occur at the zero temperature [5] . QPT be- 
ing purely quantum phenomena driven by quantum fluc- 
tuations, is associated with levels crossing and imply the 
lost analyticity in the energy spectrum. In the parame- 
ter space the points of non-analyticity, being referred to 
as critical points, define the QPT [5]. For the Hermitian 
Hamiltonian coalescence of eigenvalues results in differ- 
ent eigenvectors, and related degeneracy referred to as 
'conical intersection' is known also as 'diabolic point' 
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Since QPT is accomplished by changing some pa- 
rameter in the Hamiltonian of the system, but not 
the temperature, its description in the standard frame- 
work of the Landau-Ginzburg theory of phase transi- 
tions failed, and identification of an order parameter is 
still an open problem [8] . In this connection, an issue of 
a great interest is recently established relationship be- 
tween geometric phases and quantum phase transitions 
[9] [lOl [TTJ [12] . This relation is expected since the geo- 
metric phase associated with the energy levels crossings 
has a peculiar behavior near the degeneracy point. It 
is supposed that the geometric phase, being a measure 
of the curvature of the Hilbert space, is able to capture 
drastic changes in the properties of the ground states in 
presence of QPT [TU1 QH E3 EU 

In the context of the Berry phase the diabolic point 
is associated with 'fictitious magnetic monopole' as fol- 
lows. Assume that for adiabatic driving quantum sys- 
tem two energy levels may cross. Then the energy sur- 
faces form the sheets of a double cone, and its apex is 
called a "diabolic point" [14]. Since for generic Her- 
mitian Hamiltonian the codimension of the diabolic 
point is three, it can be characterized by three pa- 
rameters R = (X, Y, Z). The eigenstates |n, R) give 
rise to the Berry's connection defined by A„(R) = 
i(n, R|VR.|n., R), and the curvature B„ = Vr x A„ 
associated with A„ is the field strength of 'magnetic' 
monopole located at the diabolic point [T5] [TB]. The 
Berry phase 7„ = § c A„ • d"R is interpreted as a holon- 
omy associated with the parallel transport along a cir- 
cuit c pi]. 

Geometric phases and quantum phase transitions. 
- Consider the diagonalizable Hamiltonian H(X) = 
Sj=i Ei(\)\ilJi(\)){ipi(\)\, depending on the parame- 
ters A a , a — 1,2, ... ,r. Its ground state is given by 
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IV'g(A)) = <8>i=ilV'i(A)), and employing the standard for- 
mula for computing of the Berry phase, we obtain 



8 N 
M (4(A)| — |^(A)}dA a = ]T 7l (1) 

i— 1 



where 7, = § c dA^ is the geometric phase associated 
with the eigenvector |^(A)). Then applying the Stokes 
theorem we obtain 



N 



N 



^i\V a H\^ m )^ m \V b H\^i)d\ a A d\ b 



{E„ 



Ei 



It follows herefrom that the curvature FW = gL4.W di- 
verges at the degeneracy points, where the energy levels, 
say E n and E n +\, are crossing, and reaches its maxi- 
mum values at the avoided level crossing points. Thus, 
the critical behavior of the system is reflected in the ge- 
ometry of the Hilbert space through the geometric phase 
of the ground state. 

Since in the vicinity of the level crossing point only 
the two-dimensional Jordan block related to the level 
crossing makes the most considerable contribution to 
the quantum evolution, the ./V-dimensional problem can 
be described by the effective two-dimensional Hamil- 
tonian which can be obtained as follows. Let A c be 
a crossover point at which the energies E n (X c ) and 
E n +i{X c ) coalescenc. In the two-dimensional subspace 
corresponding to E n (X c ) and E n +i(X c )> we choose an 
orthonormal basis {|0),|1)} and complement it to the 
complete basis of the iV-dimensional Hilbert space by 
adding the eigenvectors IV'fc(Ac)) (k ^ n,n+ 1). 

Now, an arbitrary state \tp(t)} can be expanded as 

^{t)) = <*(t)\0) +/?(i)|l> +ES fc U,n+l^(*)IV'fc(Ac)). 

Inserting this expansion into the time-dependent 
Schrodingcr equation, we obtain the coefficients a and 
[3 as the solution of the two-dimensional Schrodinger 
equation 



;-K*)>=w e/ (A)K*)>, 



(2) 



where 



n ef (A) 



A + Z X-iY 
X + iY X -Z 



(3) 



and \u(t)) = ^ * 
are determined by 



The matrix elements in Eq. ([3]) 



A = 


~«0|F(A)|0> + 


(1|^(A)|1)), 


(4) 


X = 


~«0|F(A)|1> + 


(1|#(A)|0», 


(5) 


Y = 


5«0|ff(A)|l> - 


(l\H(X)\0)), 


(6) 


Z = 


|«0|^(A)jO> - 


(l|ff(A)|l». 


(7) 



Thus, in the neighborhood of diabolic point only 
terms related to the invariant subspace formed by 
the corresponding two-dimensional Jordan block make 
substantial contributions. The iV-dimensional prob- 
lem becomes effectively two-dimensional, and the quan- 
tum system can be described by the effective two- 
dimensional Hamiltonian H e f = XqTL + R ■ er, where 
R(A) = (X,Y,Z) (for details see [H[T9l[20]). 

The geometric phase in neighborhood of the diabolic 
point can be written as follows 
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R • dS 
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where integration is performed over the surface Ec5 2 . 
The behavior of the geometric phase described by the 
first term is independent of a peculiarities of quantum- 
mechanical system. One can consider the Bloch sphere 
as an universal parameter space for description of QPT 
in the vicinity of the critical point [21] . 

Following [TU] , we define the overall geometric phase 
of the ground state as 7 ff = (1/N) 7i- I n the ther- 
modynamical limit 7 g = J r y(x)d/j.(x), where d/i(x) is 
the suitable measure. As has been shown by Zhu [TT] on 
example of XY spin chain, the overall geometric phase 
associated with the ground state exhibits universality, 
or scaling behavior in the vicinity of the critical point. 
In addition, the geometric phase allows to detect the 
critical point in the parameter space of the Hamiltonian 
0E3Q31QI1Q3]- These works indicate that the overall 
geometric phase j g can be considered as the universal 
order parameter for description of QPT. 

The model. - The multielectron ion with in the 
crystal field has the energies of terms for d n configu- 
rations determined numerically by the Tanabe-Sugano 
diagrams [55] as a solution of the eigenvalue problem. 
Simple analytical calculations of the low energy terms 
with different spin value that is sufficient to study spin 
crossover has been done recently [23]. The crystal field 
parameter A increase linearly with pressure P. Thus 
the multielectron energies for spin Si and S2 ( E\ and 
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E2 ) are also linear functions of P. To distinguish two 
different spin states in the lattice we introduce the Ising 
pseudospin states \i) and | —i) for \d%,S\) and \df,Sl), 
where i runs over all sites in the lattice. Thus we ne- 
glect the spin degeneracy of the d™ terms but capture 
the possibility of energy level crossing that is the essen- 
tial part of the spin crossover. Then, in the basis \ + i), 
I — i) , the Hamiltonian of the system can be written as 
follows 




(9) 



where A = {E{ + E l 2 )/2, e t = {E\ - E|)/2, and 1, 
a z are the identity and Pauli matrices, respectively; the 
Hamiltonian of interaction between the spins being . 
The main contribution to the is given by the Heiscn- 
berg ecxange interaction. We consider spin crossover far 
from the thermodynamic phase transition in the para- 
magnetic phase, it allows us to simplify this interaction 
and substitute it with the effective mean field. This 
mean field is spatially uniform for the ferromagnetic in- 
sulator or two-sublattice for the antiferromagnet one. 
Examples given above [U [2j El [4] correspond to the anti- 
or ferrimanetics. In any case this mean field just renor- 
malizes the interionic multielectron energies E\ and E2 ■ 
Nevertheless this mean field results in the collective be- 
haviour of spin system in the crystal, that is why we can 
treat all spins in the space uniform states. Thus spin 
crossover at zero temperature is the transition from one 
spin ordered state to another spin ordered state. An- 
other interaction that is smaller then the exchange one 
is given by relativistic anisotropy contribution to the 
Hij . For example a spin-orbital interaction can mix dif- 
ferent spin states inside single ion, and it occurs to be 
important in our problem. 

In what follows we will consider the simplified spa- 
tially uniform model described by the following Hamil- 
tonian 



Fig. 1. Energy levels crossing 
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H = y^(A l + ea z + Xa + + A*<x_). 



where | + 1) and | - 
& z \±l ) =±\± 1 
E = \J e 2 + p 2 , we denote p 
we obtain H = 

Hi 



1) are eigenstates of the operator a z : 

u = + e/E, v 



y/l-e/E and 
|A|. After diagonalization 



^2 Hi , where 

= £+\<Pl){<Pl\ + 



r_|v3 2 )(V2| 



(13) 



and e± — Xq ± E. Due to perturbation there is a fi- 
nite gap 2p between eigenstates at the crossover point 
e = (see FigfT J. At A -> 0, u -> yT + e/\e\ and 
v — > yl — e/\e\ . When P < P c we have u — > \/2, 
v -» and \(p ± ) -> 1 + 1), \<p 2 ) -> |-1). After 
crossover (P > P c ) u — > 0, v — > y2 and \<pi) — > | — 1), 
\ifi2) — > 1 + 1)- For A ^ we can ascribe the definite 
spin to the ground term \ip2) only asymptotically. In 
order to study the geometric phase in this system, we in- 
troduce a new Hamiltonian H(V, A, ip) = U((p)WA ] '(<p), 
where U(<p) = e 1 ^" 7 * and tp is slowly varying parame- 
ter, < ip < 2tt [21] . The transformed Hamiltonian Hi 
takes the form 



n, 




(14) 



(10) 



The energy spectrum is given by e± — Ao ± y 'e 2 + p 2 , 
and the ground state energy is £_. The instantaneous 
eigenvectors are found to be 



The spin-orbit coupling A mixes the different spin states, 
and it plays the role of quantum fluctuations in our Ising 
pseudospin basis. Both Ao and e are pressure dependent. 
Further we assume e(P) = £0 — aP- The crossover takes 
place when e(P c ) — at P = P c = e^/a. 

The Hamiltonian p0[) is diagonalized by the unitary 
transformation 



sin : 



H + i) 



1)) 



l^2> = -(~V\ + 1)+U\-1)) 



(11) 

(12) 



cos • 



(15) 



where cos 8 = e 



Coupling of eigenvalues e+ 



and £- occurs at the diabolic point located at the origin 
of coordinates. 

The connection one-form associated with the ground 
state is given by A = (v,-\d\uJ) = |(1 — cos6)dip, and 
computation of geometric phase yields 7 = § c A, where 
integration is performed over the contour C on the two- 
dimensional sphere S 2 . Let us assume that the contour 
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C of integration is chosen as 8 = const. Then the geo- 
metric phase related to the ground state is 

7 = 7r(l-coB0)=7r(l--^== j. (16) 

The lost of analyticity occurs at the diabolic point lo- 
cated at the origin of the parameter space (3?A, 3A,e). 
In vicinity of the diabolic point the geometric phase be- 
haves as a step function 

= f 0, for p = 0, e -> +0 (9 -> 0) 
7 \ 2tt, for p = 0, e -> -0 (0 -> tt) 

The geometric phase 7 and its derivative d^/de ver- 
sus p, £ are plotted in Fig. O As can be observed, the 
geometric phase behaves as the step-function near the 
diabolic point, and at the diabolic point one has the 
discontinuity of the geometric phase with the gap of 2tt. 




Fig. 2. Geometric phase 7 (left) and its derivative 
dy/de (right) as a function of the Hamiltonian param- 
eters p and e. There is clear step-function behavior at 
the diabolic point p = e = 0. 

The overall geometric phase j g = (7r/-AQ^i7» can 
be written as j g = n(l + dE g /ds), where E g is the 
ground state energy per spin [21j . Besides, one can show 
that 7g = tt(1 + (<r 2 )), where (a z ) = (l/N)(i/)g\a M \^g) 
is the average sublattice magnetization per ion. In 
our model j g coincides with partial geometric phase 
7i = 1 — cos 9, therefore the step-like behavior of (a z ) 
due to high spin-low spin term crossover reported in [24j 
has the topological nature. 

Concluding remarks. - In the limit p — > the 
crossover becomes the QPT. The latter has a pure topo- 
logical nature and emerges as the quantum transition 
between the ground states with the distinct winding 
numbers [25, 26 . This quantum number, being defined 
by the geometric phase, is related to a winding number 
of the map S 1 — > U(l). Thus, the geometric phase can 
be considered as the topological order parameter in the 
spin crossovers phenomena. 
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